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Abstract 

Using first-principles based effective Hamiltonian and finite temperature Monte Carlo simulations 
we investigate cooperative responses, as well as microscopic mechanism for vortex switching, in zero- 
dimensional Pb(Zro.5Tio.5)03 nanoparticles under curled electric fields. We find that the generally 
accepted domain coexistence mechanism is not valid for toroid switching. Instead dipoles are 
shown to display unusual collective behaviors by forming a new vortex with perpendicular (but not 
opposite) toroid moment. The strong correlation between the new and original vortices is revealed 
to be critical for reversing toroid moment. Microscopic insight for the puzzling collective response 
is discussed. Based on our finding, we further describe a technological approach that is able to 
drastically reduce the magnitude of the curled electric field needed for vortex switching. 
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In the study of ferroelectric physics, one field that is currently attracting widespread 



Among FE nanostructures 



attention is low- dimensional ferroelectrics (FEs) 
that have been synthesized, FE nanoparticles are particularly interesting because the strong 
depolarizing fields along all three directions force dipoles in these particles to adopt a topo- 
logical vortex configuration, [6] which is a new form of dipole ordering [7,]. This ordering 
opens up a different class of cooperative phenomena and phase transitions in FEs of reduced 
dimension, which is a subiect that has generated long-standing interest for decades but is 



still poorly understood 



i[iol,[iil,Q,t3,Q. 



Another field of fundamental interest in ferroelectrics is to influence dipole behaviors by, 
e.g., applying electric fields 1^, 16 1 or pressure jl?!. Q. The interest stems from the fact 



that external stimuli often alter the balance among different interactions inside FEs \1{ 
thereby giving rise to new phenomena. In this context, curled electric field (namely, the 
electric field generated not by static charges, but by time-dependent magnetic field via 
V X Ecuri = —dB/dt) is an interesting route for modifying the behaviors of dipoles since the 
field possesses non-vanishing curl. The Ecuri fields can be generated on FE particles using 
focused laser beams and by varying the B-field component of the laser. But effects of curled 
fields in FE nanoparticles have never been studied before (to our knowledge). The knowledge 
as to how the interaction between FE dipoles and curled electric fields changes the phase 
properties remains entirely unknown. Technology wise, understanding the response of FE 
vortex to curled electric fields is of immense relevance, since use of bistable vortex states 
promises to increase the density of nonvolatile ferroelectric random access memories by five 
orders of magnitude. Gj Reversing the toroid moment, by applying curled electric fields, is 
a key process for making this possible. Recently, a study of using static electric fields to 



influence the formation of vortex in FE particles was reported[20], where the system is first 
heated to high temperatures so that the original vortex disappears, and then the paraelectric 
particle is annealed under the influence of external static fields. By adopting different 
external fields, vortices of varied orientations were obtained after annealing. Obviously this 
is not a switching process in which a constant temperature is maintained. 

Independent of both subjects given above, a third topic of importance is to understand 



switching mechanism in FEs. In bulk FEs, coexistence of 



accepted as the mechanism for polarization switching. 



and P domains has been widely 



15| More specifically, the switching 



begins by nucleation of a -P domain of reversed polarization; crystal defects (domain walls. 
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2l| Subsequently, the -P domain grows 



twin boundaries, etc.) often are the nucleation sites, 
while the original P domain shrinks under external fields, leading eventually to switching. 
In FE nanoparticles, formation of domain walls is less likely and energetically unfavorable. 
However, the vortex phase in FE particles does have an intrinsic topological defect — the 
vortex core. Based on the fact that a large strain energy exists near the center of the vortex 
where dipoles turn their directions sharply, vortex cores are thus less stable. When imposed 
under a switching curled field, dipoles near the vortex core are likely to fiip first to form a 
-G domain of reversed moment. This -G domain coexists with and gradually conquers the 
G domain, leading to switching. 

Here we perform ab initio based studies to investigate for the first time the collective re- 
sponse to curled electric fields, and the vortex switching mechanism, in ferroelectric nanopar- 
ticles made of Pb(Zro. 5X10.5)03 solid solution. The purpose of this Letter is (1) to demon- 
strate that the generally accepted mechanism of domain coexistence is invalid for vortex 
switching, and in fact, vortex reversal occurs by following an interesting Gz-Gy-{-Gz) phase 
transformation sequence; (2) to report the unusual collective phenomena and energetics in 
OD EEs in response to curled electric fields. Based on the key insight obtained in this study, 
we further propose and demonstrate an alternative approach (which could be technologically 
very important) to overcome the energy barrier of vortex switching. It should be pointed out 



that this work is profoundly different from a previous study[24] where effects of homogeneous 
electric fields in EE particles were investigated. Homogeneous electric fields, possessing no 
curl, can only transform a vortex phase into a ferroelectric phase of polarization. As a result, 
the switching of toroid moment, which is the key subject of this study, was not (and can 
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not even be) addressed in Ref. 

Technically, we use first-principle derived effective Hamiltonian 23|, |2J] and finite tem- 
perature Monte Carlo (MC) simulations. Besides internal energy U, the Hamiltonian also 
consists of a — E(rj) ■ pi term,|l^ describing the coupling between external electric field 
E(r) and local dipoles pj. Nanoparticles of disk shape are considered, with cylindrical 2;-axis 
parallel to the crystallographic [001] direction. Switching simulations under curled electric 
fields are performed at 64K where the vortex phase is stable. A fixed acceptance rate of 0.25 
and 10000 MC sweeps are used throughout all simulations. 25| We will mainly present re- 
sults for a disk with diameter d=19 and height h=14 (both quantities are measured in lattice 
constant a=4A of pseudocubic bulk), since the diameter (7.6nm) of this disk falls into the 
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size range (5-50nm) of experimentally fabricated nanoparticles 26|]. In this work we use MC 
rather than molecular dynamics (MD) to study switching dynamics, based on the following 
considerations, (i) In both methods the basic quantity governing the dynamics is the same, 
which is the energy. The difference is that the energy is reflected in MD by using forces, 
while in MC the energy itself is used in determination of next configuration, (ii) Monte 
Carlo step can be quantified as a "quantum" time, and the equivalence of two methods was 



rigorously demonstrated in Ref, 



271 . where Monte Carlo simulations and Langevin dynamics 



yield excellent agreement, (iii) Vortex reversal in FE particles is essentially a thermally 
activated process (which is slow). MD simulations will be time consuming. But computa- 



23 is extended here to include 



tions of MC simulations are moderate. The method in Ref|! 
both static homogeneous electric field and curled E^uri field, i.e., E(r) = Eh + Ecuri(r)- 
The latter is given as Ecuri = | S x r, where S = Sez measures the vorticity of the field. 
With the Eh and Ecuri fields acting together, the interaction energy — E(ri) ■ pj be- 
comes -N (Eh ■ (p) + S • G), where N is the number of 5-atom cells in nanoparticles, (p) 
the average dipole moment, G = ^ toroid moment. Vector S thus is the 

thermodynamically conjugated field for toroid moment G. 

Figure [1^ depicts the collective behaviors of the toroid moment G, evolving during sim- 
ulation as a function of MC sweeps (denoted as n, in units of 200), in the d=19 disk under 
a curled field of S=0.25 mV/A^ that has an opposite vorticity with respect to the initial FE 
vortex. One key finding in Fig{T^ is that the vortex reversal process is predicted to consist of 
three evolution phases in which the system displays distinct dipole behaviors. These phases 
are, when MC sweep n is below ni=30 (to be named as phase I), n is between ni=30 and 
n2=43 (to be referred to as phase II), and n>n2 (phase III). At the initial instant of the 
simulation, the vortex has only a toroid moment, while the G^ and Gy components are 
null. As the curled field is turned on, most notable results occurring in Fig{T^ are: (1) The 

moment remains remarkably stable and subjects to only marginal decrease throughout 
phase I. However, the system does not idly wait; another critical activity is taking place 
in phase I. That is, the FE nanoparticle is building its G^ moment which increases appre- 
ciably and becomes very significant (with a value of 15 eA^) in the end of phase I. (2) As 
the system enters phase II, dramatic difference occurs. The G^ moment starts to decline 
sharply, reverses its direction, and then is stabilized at its negative saturated value at n2. 
(3) While Gz switches its direction in phase II, the G^, moment undergoes a complex evolu- 
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tion, by first rising (with a different slope than the preceding process in phase I) and later 
dropping swiftly to null. Precisely at the moment when G2 is zero, Gy reaches its maximum 
magnitude. 

From simulations, the key toroid component during switching is G^, G^, and -G^ for the 
phases of n<ni, ni<n<n2, and n>n2, respectively. The reverse of FE vortex therefore occurs 
by undergoing an evolution sequence of Gz Gy ^ (— G^), in which appearance of the 
lateral Gy moment is a critical bridging process. Formation of this lateral vortex is unusual 
(and puzzling), for neither the initial vortex of the system nor the curl of the external 
field have nonzero components in the lateral xy directions which are capable of causing the 
observed collective behaviors. To verify the above switching sequence, calculations with 
varied initial dipole configurations or temperatures or slightly different Ecuri fields were 
performed. All confirm the occurrence of a lateral vortex as the bridging phase. We also 
find that the toroid moment of the intermediate lateral vortex can point at any direction 
within the xy plane, since x- and y-axis are equivalent. In simulations of Figil^, the lateral 
vortex happens to be mainly along the y-axis. 

To explain why the Gy toroid moment forms and how it influences the reversal of the G^ 
component, we now provide microscopic insight on the field- induced dipole behaviors. Figl2] 
shows the snapshots of local dipole distributions corresponding to selected MC sweeps. More 
specifically, Fig|2^-d describe the development of the Gy toroid moment, while Figl2^-h show 
the evolution of the inplane G^ vortex. Our discussion shall follow the sequential order of 
evolution. Initial response to the curled electric field is subtle, and begins with the dipoles 
located at the vortex core. Common wisdom tells us that, under the Ecuri field, the dipoles 
near the vortex core should flip into opposite directions to form a -G domain. Interestingly, 
this does not take place in our simulations. We find instead that these dipoles rotate 
towards the 2-axis (Fig 12^). The dipole rotation is, however, energetically less favorable, 
since it generates a z-direction depolarization field that can not be compensated by the 
curled field. Note that this situation is unlike the case when a static homogeneous electric 
field is applied in the z-axis, where the homogeneous field can overcome the existence of the 
depolarization field. The delicate balance between the internal depolarization field and the 
external curled field leads to the dipole configuration in Figj2^, where the dipoles in the lower 
part point downwards as a result of the curled field while the dipoles in the upper part point 
horizontally to reduce the depolarization field. As evolution proceeds to n=20 (Figj2b), 
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dipoles of the lower part rotate also horizontally, but opposite to the dipole direction in the 
upper part, to avoid lateral depolarizing field. Consequently, the dipole pattern developed in 
Fig 12)3 explains how the nucleation seed of lateral Gy vortex is formed (which is important). 
Another key feature of the dipole pattern in Fig|2)3 is that the left and right arms of the 
vortex are largely absent. Interestingly, as the Gy vortex is nucleated at n=20, the inplane 

vortex is nearly intact (see the G^ vortex in Figj2^, also at n=20), which is consistent 
with the result in Fig{T^ where the G^ moment remains robust throughout phase I. 

The system behaves differently when evolving into phase II, highlighted in Figl2]f by the 
dramatic disappearance of inplane dipole components in regions A and B. More precisely, 
dipoles in region A rotate toward the negative z-axis, while dipoles in region B toward the 
positive z-direction. Disappearance of the inplane components for dipoles in A and B regions 
leads to three marked effects: (i) the onset of a sharp decline of the Gz moment, as seen 
in Fig{T^ at the beginning of phase II, (ii) the core of the G^ vortex to shift away from the 
cylindrical axis (Figl2]F), and (iii) a rapid increase of the z-component magnitudes for those 
dipoles located at the left and right arms of the lateral Gj^ vortex, as revealed by contrasting 
FiglJt with Figj2b. It should be recognized that nucleation of the lateral G^^ vortex, right 
in the form of Figl2]D, is critical for the destruction of the original G^ vortex, since the 
oppositely oriented horizontal dipoles at the upper and bottom surfaces in Figl2]D produce 
only a small strain energy when the dipoles in regions A and B rotate into the z direction, 
and simultaneously, no significant increase in the depolarizing field occurs along the z-axis. 

As the Gz vortex continues to be uncurled, the G^ moment becomes zero at n=33 (see 
FigH^). A vanishing G^ moment does not mean disappearance of the xy components of all 
dipoles, however. Instead, calculation result in Figl2^ reveals that the dipole components 
on the xy plane evolve, interestingly, into a stripe pattern, which also yields zero G^. The 
stripe is formed due to the fact that the curled field continues to push the G^-vortex core 
in Fig|2]F to the (lower) surface. Starting with the stripe in Figl2^, dipoles in the upper part 
rotate and/or fiip their directions, forming the critical embryo of the reversed -Gz vortex at 
n=34 (see Figl2)i). After the nucleation of this -G^ vortex, the lateral Gy vortex starts to 
disassemble itself by rotating its dipoles at both arms back to the xy plane (see Figjjji at 
n=40), which eventually leads to the development of a full and reversed -Gz vortex. 

From the above striking microscopic insight, a novel mechanism that governs the switch- 
ing of FE vortex is thus obtained. That is, vortex reversal occurs, not by the coexistence of 
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and -G^ domains, but by formation of a new and lateral vortex. In fact, we find from the 
entire evolution process in Figl2]that, the G^ and -G2 domains never coexist at any time in 
the simulations. Only after the original G^ vortex is dismantled does the nucleation of the 
-Gz vortex start to form (FigI2)i). We further numerically find that the magnitude of the 
total toroid moment |G| = y^G^ + Gy + is almost a constant throughout the switching 
(see the diamond symbols in Fig{T^), demonstrating that vortex reversal occurs by obeying 
the conservation of the |G| moment. 

We now study the energetics of vortex switching. Here we are interested in internal 
energy U, rather than free energy F=U-^^ E(rj)-pj, since the former yields important 
knowledge concerning the energy barrier. The internal energy is given in Fig{T^, revealing 
new observations that do not appear in collective behaviors. First, the highest energy barrier 
occurs around n=30, which corresponds to the instant when the Gy vortex has formed and 
the inplane G^ vortex just starts to annihilate its xy components. Height of the energy barrier 
is 2.4 meV per 5-atom cell. Second, and interestingly, a local energy valley appears at n=33, 
showing that the fully developed Gj^ vortex (with G2=0) is actually a stable structure even 
if the curled field is removed. Third, the two energy barriers surrounding this energy valley 
are asymmetric, with a lower barrier of 0.8 meV at n=35 (favoring the formation of the 
-Gz vortex). To examine how the energy barrier may vary with temperature and size, we 
have performed switching calculations at 160K (as compared to 64K in Fig{T^), for the same 
d=19 and h=14 dot. The energy barrier is found to be 2.1 meV at 160K, close to the value 
of 2.4 meV at 64K. This is consistent with the fact that the energy barrier is determined 
mainly by depolarization field, not by temperature. We have also performed calculations 
for other d diameters, all showing that formation of the lateral vortex is a key step for 
switching the G^ vortex. We further find that depending on the h/d ratio, the number of 
lateral vortices — that occur at the switching moment (i.e., when G^ equals zero) — differs. 
In the d=9 and /i=14 disk with h/d ^2, two Gy vortices and one G^^^ vortex (with dipole 
configuration similar to the structural B phase in Refill) take place at the switching moment. 
The energy barrier for the d=9 disk is found to be 1.7 meV at 64K. 

Interestingly, we previously found that homogeneous Eh = Ehe^ electric fields can also 
give rise to a lateral vortex |22[|, despite that homogeneous fields and curled fields are drasti- 
cally different in nature. Based on the key findings that (1) the and E^uri fields both are 
capable of inducing the Gj^ vortex, and (2) formation of the Gy vortex is the major energy 
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barrier for toroid switching (as demonstrated in this study), we further propose an alterna- 
tive approach to switch FE vortex, by means of a combined action of the Eh and Ecuri fields. 
More precisely, an Eh field of sufficient magnitude is applied to precipitate the formation of 
Gy vortex, and an Ecuri field is then added on (while the Eh field remains) to generate the 
reversed vortex. A significant advantage of this approach is, while a strong curled field 
is hard to achieve since it requires fast oscillation of alternating magnetic fields, a strong 
Eh field of lO^V/m order can be easily generated on nanometer scale. We demonstrate the 
approach using the same d=19 disk, by placing it first under a homogeneous Eh=1.9V/nm 
field (which quickly induces a vortex), and then under a curled Ecuri=0.04mV/A^ field 
for ten thousand MC sweeps. Fig{T]D describes the evolution of the toroid moment after the 
curled field is turned on. Drastic differences are evident in FiglDo as compared to FigH^: (1) 
When the curled field is switched on, the Gj, value in Fig{TjD is already notably large. This 
Gy moment is formed by the homogeneous field, confirming the result of Ref]22l. (2) The 
long period of phase I — that lasts five thousand MC sweeps in FigilK — no longer exists in 
FigHb. Instead, G^ starts to decline sharply in FigHb as soon as the curled field is imposed. 
(3) The magnitude of the curled field used in FigHb is considerably reduced, by a factor of 
600% as compared in Fig{T^. This reduction factor can be further increased when a larger 
Eh field is used. Analysis of the internal energy (solid symbols in Fig{T]D) shows that the 
energy barrier almost disappears, and the curled field acts to guide the formation of the 
reversed vortex. 

In summary, response of FE particles to curled electric fields was studied for the first 
time. Simulations revealed a novel mechanism that governs the switching of FE toroid 
vortex. A curled electric field, of only z-axis vorticity momentum, was shown to lead to the 
unexpected formation of a lateral vortex, manifesting macroscopically the balance of various 
microscopic interactions. We demonstrated that the strong correlation between the original 
vortex and the new lateral vortex plays a critical role for vortex reversal. On the other 
hand, the -G and G domains never coexist during switching, and the switching mechanism 
by coexisting domains is not valid. Moreover, based on the findings that formation of lateral 
vortex is the major energy barrier for vortex switching, and that this lateral vortex can also 
occur by use of homogeneous field, we described an effective approach that is able to reduce 
substantially (by 600%) the magnitude of the switching curled field. Finally, we believe that 
the microscopic insight and energetics (which are difficult to obtain in experiments) will also 
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be 



of immense value. 
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10 20 30 40 50 
Monte Carlo sweeps (x200) 

FIG. 1: Evolution of toroid moment G and its magnitude |G| (empty symbols, using the left 
vertical axis) and internal U energy per 5-atom cell (solid dots, using the right vertical axis) in a 
d=19 nanodisk: (a) under a S=0.25mV/A^ curled electric field; (b) under the combined action of 
a Eh=1.9V/nm homogeneous field and a S=0.04mV/A^ curled field. For clarity of display, is 
plotted after multiplying -1. 
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FIG. 2: Dipole configurations in a d=19 nanodisk under a curled field of S=0.25mV/A^ at different 
MC sweeps. Upper panel (a-d): dipoles on the central xz cross section revealing the evolution of 
the Gy vortex, at the following MC sweeps (a) n=10, (b) n=20, (c) n=30, (d) n=40. Lower panel 
(e-h): dipoles on the central xy cross section displaying the evolution of the vortex, at (e) n=20, 
(f) n=30, (g) n=33, (h) n=34. All (n) numbers are in units of 200 MC sweeps. Arrows show the 
magnitudes and directions of local dipoles, projected on the corresponding plane. For convenience, 
(b) and (e) corresponding to the same n=20 sweep are given on the same column [so are (c) and 

(f)]- 
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